Aharonov-Bohm ring with fluctuating flux 
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We consider a non-interacting system of electrons on a clean one-channel Aharonov-Bohm ring 
which is threaded by a fluctuating magnetic flux. The flux derives from a Caldeira-Leggett bath of 
harmonic oscillators. We address the influence of the bath on the following properties: one- and 
two-particle Green's functions, dephasing, persistent current and visibility of the Aharonov-Bohm 
effect in cotunneling transport through the ring. For the bath spectra considered here (including 
Nyquist noise of an external coil), we find no dephasing in the linear transport regime at zero 
temperature. 

PACS numbers: 73.23.-b, 73.23.Hk, 73.23.Ra, 03.65.Yz 
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I. INTRODUCTION 

In the present work, we consider a simple theoretical 
model system of non-interacting spinless electrons which 
are restricted to move in one dimension around the cir- 
cumference of a clean, one-channel Aharonov-Bohm ring. 
The ring is threaded by a magnetic flux which fluctu- 
ates around some average value (see Fig. This may 
lead to dephasing of the electron motion on the ring, 
apart from other effects like renormalization of the elec- 
tron masses and introduction of an effective coupling 
between the electrons. We treat the full dynamics of 
the fluctuating flux coupled to the electron system in a 
selfconsistent manner, rather than prescribing an exter- 
nal stochastic time-dependent classical field. In order to 
achieve this, the flux is taken to be the sum of the normal 
coordinatespof a Caldeira-Leggett type bath of harmonic 
oscillatorstlcl. The fluctuations couple to the electrons via 
the vector-potential term in the kinetic energy. As an im- 
portant special case for the bath spectrum we treat the 
Nyquist noise that may be due to the equilibrium current 
fluctuations in the external coil producing the flux. 

Equilibrium and transport properties of this model sys- 
tem are analyzed for the cases of zero and finite temper- 
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FIG. 1: The model situation: A fluctuating flux leads, via a 
time-dependent vector potential, to a fluctuating force for the 
electrons on the Aharonov-Bohm ring 



atures, taking into account the coupling to the bath and 
the Pauli principle with respect to the electron system. In 
particular, we discuss the single-particle and two-particle 
Green's functions, level widths, energy shifts and dephas- 
ing times, and the reduction of the persistent current due 
to the fluctuations. In each case, the dependence on the 
coupling strength between system and bath and on the 
low-frequency spectral properties of the bath is exam- 
ined. Aharonov-Bohm interference observed in cotunnel- 
ing through the ring is discussed in order to analyze the 
coherence properties of the electrons on the ring under 
the influence of the fluctuating flux in a transport situ- 
ation. As a result of our calculation, we find that the 
fluctuations do not lead to dephasing in the linear trans- 
port regime at zero temperature. 

The single-particle version of this model has been con- 
sidered before in Ref. ^ in order to determine whether 
persistent currents in a normal metal ring may be de- 
stroyed by coupling to an Ohmic bath. Whereas the 
authors oL Ref. ^ used the Feynman- Vernon influence 
functionalu, we will apply a different, more direct ap- 
proach. We emphasize that our analysis is restricted to 
baths weaker than the Ohmic bath (see the discussion af- 
ter Eq. (|5^)). The possibility of spontaneous persistent 
currents was investigated (and ruled out) in Ref. ^ using 
a Luttinger liquid picture for the electrons and taking 
into account their electromagnetic self-interaction. De- 
phasing of a single electron going around the two arms 
of an Aharonov-Bohm ring has been considered both in 
Ref. 1^, using the influence functional, and in Ref. ^, using 
a semiclassical picture. In the latter article, the connec- 
tion between phase fluctuations and the trace left by the 
system in the environment was emphasized. More re- 
cently, the question of dephasing in mesoscopic systems 
has received renewed attention due to a set of weak- 
localization measurements which have shown a-|Satura- 
tion of the dephasing time at low temperaturesElB. Mo- 
tivated in part by these puzzling flndings, the authors of 
Refs. p^yil]Jl^ considered a ring containing a single quan- 
tum dot with fluctuating gate voltage and obtained the 
properties of the quantum-mechanical ground state (in 
particular the persistent current). A strong influence of 
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external nonequilibrium noise on the persistent current 
in a disordered quasi-lD ring has been found recently 
in Ref. The effects of a phase-breaking scatterer in 
the many-particle situation, where the Pauli principle be- 
comes important, have been discussed in Rcf. Very 
recently, dephasing in a mesoscopic Mach-Zehnder type 
interference setup has been analyzed in Ref. [T5[ 

In the following section, we will define the model and 
discuss some simple consequences as well as some features 
which cannot be included in this system easily. Then we 
give a short qualitative discussion of dephasing for the 
simplified case of a classical fiuctuating fiux (represented 
by a random process) . Similar considerations are applied 
to the calculation of the Green's function for a single elec- 
tron on the ring, both with classical and quantum fiuc- 
tuations. The resulting energy shifts and level shapes 
are analyzed in some detail, since these results can be 
taken over to the many-particle calculation of the single- 
particle and two- par ticle Green's functions which is pre- 
sented in section |^ together with the evaluation of the 
grand canonical partition sum and the persistent current. 
After discussing the physical meaning of the dephasing 
produced by the Nyquist noise at low temperatures, we 
turn to an analysis of Aharonov-Bohm interference in a 
cotunneling transport measurement. 



II. THE MODEL 

The Hamiltonian of the system of electrons on the ring 
is given by 



P~ 9<t> 
2m 



-4'p + Hbath , 



(1) 



where Hbath is the Hamiltonian of the set of uncoupled 
oscillators representing the bath: 



H, 



bath 




(2) 



The possible values of the electron momentum p are 
quantized due to the finite circumference L of the ring: 
p — 2Tm/L with an integer n. (Note that here and in the 
following, we will put h = ks = I) The term gcf) in the 
kinetic energy of the electrons is due to the coupling to 
the vector potential which is proportional to the fiuctu- 
ating flux, (j) represents this flux (up to a constant factor) 
and is assumed to be given by the sum over the oscillator 
normal coordinates: 



(3) 



The prefactor in this definition has been chosen such 

that the autocorrelation function ^0(i)(^(O)^ of (f) has a 

well-defined limit if the number Nose of oscillators tends 
to infinity while the spacing of frequencies tends to zero 
like 1/Nosc- This is the "thermodynamic limit" of an infi- 
nite bath which is necessary to describe truly irreversible, 
dissipative dynamics. The quantity g is the coupling 
strength between bath and electrons. It incorporates the 
electron charge and the circumference of the ring, since 
the line-integral of the vector potential around the ring 
gives the flux. Any external static magnetic fiux $ has 
to be added in the kinetic energy expression. 

We assume the interaction between system and bath to 
be sufficiently weak, such that the bath may be treated as 
linear in a good approximation (like it is usually done in 
the theory of quantum dissipative systems, see Ref. for 
a more detailed discussion). Apart from this assumption, 
the expression Eq. (||) used for the fiuctuating flux is still 
completely general. We are free to choose the frequencies 
of the bath oscillators to obtain any desired correlation 
function of 0, which is the only quantity that affects the 
dissipative system dynamics in this model. Note as well 
that the coupling of the velocity to-a vector potential as- 
sumed in Eq. (y) has been showr£3 to be equivalent to 
a Caldeira-Leggett model with the usual coordinate-type 
coupling for a particle moving on an infinite line. Phys- 
ically, the coupling used here is the natural choice for a 
situation in which the vector potential is linearly related 
to the fluctuating current in an external coil. The fluctu- 
ations (and linear response) of the current then determine 
the bath correlator discussed in the next section. 

The correlator of will determine the dephasing rate 
and other important quantities via its low-frequency 
properties. For the discrete set of oscillators and for the 
continuum limit, respectively, it is given by the following 
expressions, where the average (•) is taken with respect 
to the unperturbed set of oscillators: 



N, 



mm 



coth cos (ujjt) — i sin (ujjt) 



dLuC{uj) 



coth (^^) cos {ujt) — ism {tot) . (4) 



This deflnes the spectral function C{u!) which we will 
use to characterize the bath spectrum. In terms of the 
discrete set of frequencies, it is given by C{uj) = iVj^J. • 
J2j S {u! — U!j) {2Mujj)~^. Note that the coth is equal to 
2n{Lu) + 1, where n(uj) is the Bose-Einstein distribution 
function. The special case of Nyquist noise is obtained by 
the requirement that for high temperatures, the spectrum 
of fluctuations of the magnetic flux (i.e. of 0) is white. 
Since coth {ijj/2T) = 2T/uj for T > w, this means C{ijj) oc 
Lo (for small uj). 

Before we proceed to the calculations, we will point 
out some simplifying features of this situation as well 
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as some important aspects of the dephasing problem in 
degenerate Fermion systems that are beyond the scope 
of this model. 

The magnetic flux is assumed to thread the ring in such 
a way that the situation is axially symmetric with respect 
to the axis which goes through the center of the ring and 
is perpendicular to its plane. In this case, we can choose 
the gauge such that the vector potential is everywhere 
tangential to the ring and of constant magnitude around 
the whole circumference. The same holds for the electric 
field, which is given by the time-derivative of the vec- 
tor potential. This is analogous to the Caldeira-Leggett 
treatment of one-dimensional quantum Brownian motion 
of a free particle, with the formal difference that in our 
case the force is derispd from a vector potential instead 
of a scalar potentialBE^. It is the choice appropriate for 
a system with periodic boundary conditions, where the 
quantization of momenta, the Aharonov-Bohm effect and 
persistent currents play a role. Note that under different 
circumstances the assumption of a force which is constant 
in space is only valid within the dipole approximation for 
a particle that is restricted to move in a well- localized re- 
gion of space, because otherwise the finite wavelength of 
the bath modes (phonons etc.) becomes important. 

Since the coupling between system and bath is via the 
momentum, which commutes with the electron Hamilto- 
nian and therefore is a constant of the motion for the orig- 
inal uncoupled electron system ( "diagonal coupling" ) , 
some features of this model are very simple: In spite of 
the interaction with the bath, the momentum of a parti- 
cle will stay constant, only its velocity and kinetic energy 
can fluctuate. This simplifies the many-electron prob- 
lem as well: Although the fluctuating force influences 
the center-of-mass motion of the electrons and introduces 
some kind of "effective interaction" between them, the 
occupation of different p-states cannot be changed by the 
bath. Note that this simplification would be spoiled if 
one takes into account impurities and/or a coupling that 
depends on the position. For example, the latter would 
arise if one considered an arbitrary fluctuating electro- 
magnetic field or the electric field between the plates of 
a capacitor, which is a constant vector field in space but 
is not constant with respect to its projection onto the 
direction of motion of the electrons on the ring. Other 
situations where the coupling depends on position include 
interaction with phonons or a localized spin on the ring. 

In its single-particle version, the Hamiltonian given 
above also arises in the discussion of dephasing for a 
charged island, if the coupling to the bath (e.g. a fluc- 
tuating gate voltage) is purely diagonal in the system's 
eigenbasis and can only lead to fluctuations in the energy 
levels. Note, however, that questions of interest here like 
tunneling into the ring or features of the many-particle 
system have no natural counterpart in that rather simple 
situation. 

Although in our model all electrons are coupled to the 
same flux, which introduces a kind of effective interaction 
between them, the decay rates of Green's functions will 



not show any dependence on the distance to the Fermi 
surface, in contrast to the usual behaviour of interacting 
Fermi systems. This is due to the diagonal coupling be- 
tween system and bath, which means that there are no 
energy-relaxation processes which change the occupation 
numbers of the electrons and which would feel the restric- 
tion by the Pauli principle. A related question arises in 
the study of dephasing in degenerate Fermion systems: 
If the coupling is not diagonal in the electrons' (single- 
particle) eigenstates, the system variable which couples 
to the bath (in our case the momentum) carries out fluc- 
tuations itself. In a semiclassical single-particle calcula- 
tion, these fluctuations pick oiijt the high-frequency com- 
ponents of the bath spectrumEi. Therefore, according to 
such a calculation, there is dephasing even in the case 
of a bath spectrum that vanishes at low frequencies and 
this implies that at low temperatures, the high-frequency 
zero point fluctuations of the bath contribute heavily to 
dephasing in this picture. However, if the electron sys- 
tem is nearly degenerate, many of its transitions will be 
blocked by the Pauli principle so that such an effect will 
be strongly suppressed. Although we cannot investigate 
this point in our model, related considerations will occur 
in our discussion of the cotunneling setup in section M. 



III. SINGLE-PARTICLE PROBLEM 

In this section we will discuss the problem of a single 
particle on the ring, both semiclassically and quantum- 
mechanically. The results can support the understanding 
of the next section, which is devoted to the many-particle 
situation. 



A. Semiclassical analysis 

Consider two wave packets traversing the left and right 
arm of the Aharonov-Bohm ring with constant velocity 
V and meeting again after some time t = L/ {2v) at the 
opposite end. The resulting interference pattern depends 
on the total phase difference between the two paths. In a 
semiclassical calculation, the phase difference is produced 
by the vA-term in the Lagrangian of the particle, and it 
is given by 

ip = 2- f vA(t')dt'. (5) 
c Jo 

The factor 2 arises because the phases are equal up 
to a change in sign. A{t') gives the time-dependence of 
the fluctuating vector potential which is assumed to be a 
classical Gaussian random process with zero mean in the 
high-temperature limit considered here. The visibility 
of the interference pattern will be suppressed due to the 
fluctuations in the phase ip. Since ip is a Gaussian random 
variable, we obtain for the suppression factor 
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(6) 



In our model eAjc is equal to g(^, where we will treat 
as a classical fluctuating field by taking into account 
only the real (symmetric) part of the correlator (0) in the 
high-temperature limit. Then the variance (^^j) of the 
phase becomes: 



dti / dt2 {<j){h)(j){t2)) 



(7) 



If the traversal time t is much larger than the corre- 
lation time of the fluctuations in cj), we may apply the 
following standard approximation: 



Here and in the following quantum-mechanical calcula- 
tion, we will neglect this term, which does not couple to 
the momentum (cf. discussion in section IV E). With 
this approximation, the Green's function is given by an 
expression which again involves the correlation function 
of 0: 



iGfit) « exp 



-I 1 

2m 



X exp 



1 (9P 

2 \m 



dh / dt2 mimh)) 



(10) 



At this point, the discussion given above applies. In 
particular, for Nyquist noise and finite temperatures, the 
Green's function decays exponentially with a rate pro- 
portional to g^p^T. 



+ 00 



(^2) « t ■ Ag\^ / dt' {mm) 



t ■ Ag'^v'^ ■ AttT 



(8) 



This means that in the case of Nyquist noise (C oc lo) 
we obtain a finite "dephasing rate" that grows linearly 
with temperature. Note however, that the time t intro- 
duced here cannot grow without bounds but is fixed by 
the circumference L and the velocity v. This calcula- 
tion already shows that the dephasing rate will vanish 
for u ^ or T — > 0. For a bath that is weaker at 
low frequencies (C cx with a > 1) we do not ob- 
tain a suppression factor that decays exponentially with 
time, hence the dephasing rate is always zero. We defer 
a more detailed discussion of the various cases to the full 
quantum-mechanical treatment below. 

The physical interpretation of this result is clear: The 
fluctuating electric field oc A leads to a fluctuating ve- 
locity (X A, so that the random shift in the interference 
pattern is Ax oc / A{t') dt' . The interference pattern will 
be completely washed out once the spread in Ax becomes 
comparable to the wavelength A oc 1/v. This coincides 
with the criterion {f"^) ~ 1. 

Very similar considerations arise in the calculation of 
the Green's function of a single electron on the ring. The 
retarded Green's function may be approximated semi- 
classically by averaging the amplitude for propagation of 
the electron under the influence of the fluctuating flux: 



0(<)({*p(t),*t(O)} 



exp 



2m 



dt' 



(9) 



The exponential contains a quadratic term {g4>)'^ which 
does not represent a Gaussian random variable, so that 
formula (^ cannot be applied to perform the averaging. 



B. Single-particle Green's function: Quantum case 

In the following, we will calculate and discuss the 
quantum-mechanical Green's function of a single electron 
on the ring. The extension to the many-particle case will 
be given in the next section. The single-particle density 
of states (DOS), which is given by the imaginary part of 
the Fourier transform of , is a measurable quantity, 
as it can be revealed by tunneling into the ring. 

We imagine a situation where the ring is empty and a 
single electron is inserted, so that the retarded Green's 
function is 



iGi'it) = e{t) = e{t) 



iHt 



(11) 

The average is a thermal expectation value with re- 
spect to the unperturbed bath of oscillators, correspond- 
ing to the situation without any particle on the ring. 
Obviously the bath cannot change the occupation of the 
different momentum states. Therefore, we only have to 
take into account that the time evolution between cre- 
ation and destruction of the electron is governed by the 
Hamiltonian which contains the kinetic energy term that 
couples to the bath via the flux (p. Through this cou- 
pling, the introduction of the particle into the ring per- 
turbs the bath oscillators. In anticipation of the many- 
electron case, we will denote by if [{pj }] the Hamiltonian 
for a flxed number of occupied momentum states, given 
by the set {pj}- This operator only acts on the bath. 
Then the matrix element of the time-evolution operator 
with respect to the electron system is given by the fol- 
lowing expression, for a Slater determinant belonging to 
the configuration {pj}'. 



(fe}|< 



iHt 



\{pj}) = e-'-f^[feM*, 



(12) 



In particular, without any electrons we have -ff[0] = 
Hbath- For the Green's function considered here, this 
leads to: 
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= ( T cxp 



2m m Jq 



13) 



This expectation value can be interpreted as the (ther- 
mally averaged) overlap between the initial bath state 
evolved once with and once without presence of a parti- 
cle on the ring. In the second line, we have introduced 
the interaction picture with respect to H^ath and dropped 
the term quadratic in ^0 from the kinetic energy (com- 
pare the discussion above). Since </> is a bosonic vari- 
able (i.e. linear in the oscillator normal coordinates), 
Wick's theorem can be applied to the evaluation of this 
time-ordered thermal average, using a linked-cluster ex- 
pansion. It leads to an expression completely analogous 
to the one used above for the classical Gaussian ran- 
dom process, see (^) and (10). The difference consists 
in the replacement of the classical correlator by the ther- 
mal time-ordered expectation value: 



IK (j){t') 



exp 



T exp 



(14) 



At present, k = gp/m, but the same formula will be 
used below with other values for n. In contrast to the 
classical correlator, (T(l){ti)(j){t2)^ is complex and will 

lead to an energy shift in addition to a decay of the 
Green's function. Using 



duo C{uo) [(2n(w) + 1) cos {uot) ~ i sin {lo \t\)] (15) 



the presence of the particle, so the overlap between initial 
and final bath state is diminished and the Green's func- 
tion decays. The first integral produces a negative energy 
shift, corresponding to the formation of a new interacting 
state of particle and bath. Since this shift is proportional 
to , it is equivalent to an enhanced effective mass m* . 

Combining the energy-shift with the initial kinetic en- 
ergy p'^ /2m, we obtain the following expression for m*: 



1 If g 
= — 1 - 2— 

m* m \ m 



2 , C{cj) 

dio 

^ 



m 



(1-0- (17) 



Alternatively, the effective mass can be calculated from 
the initial Hamiltonian (for a single particle) by us- 
ing the term gpcji/m from the particle's kinetic energy 
to introduce a momentum-dependent shift into the os- 
cillator potential energies of Eq. (^): Qj ^^ Qj ~ 
gp{mM(jjj \/Nosc)~^- The resulting term quadratic in 
p must be compensated for, which yields the change in 
mass: 



1 

m* 




9 



mM N, 



-y- 

OSC "^o 



(18) 



Since we have neglected the term (^5^^ in the preced- 
ing derivation, the effective mass displayed here is only 
correct in lowest order with respect to g^. A full calcu- 
lation yields to* = to (1 -|- ^). 

The real part of the second integral in Eq. (^6|) gives 
a negative contribution, corresponding to a suppression 
in magnitude of the Green's function. If the bath is rel- 
atively weak at low frequencies (C oc with a > 3), 
the long-time behaviour is simple: The decay saturates 
at i ^ oo. 

In this case, the Fourier transform of the Green's func- 
tion. 



to evaluate the double time-integral, we obtain for the 
exponent of Eq. (n4) : 



2 / , C(lo) 







P diu [n {\lu\) + d{uj)) (e-^"* - l) . (16) 

' —oo 



The step-function d{Lo) corresponds to the zero-point 
fluctuations. In the remainder of this section, we will dis- 
cuss the behaviour of the Green's function derived from 
( p^ ) at short and long times, for different bath spectra. 

At short times, the integral in the second line compen- 
sates the first integral, so that the exponent begins to 
grow like instead of t. Physically, this means that at 
short times the particle has not yet influenced the bath 
and its energy is still given by the bare energy p^ /2m. 
At later times, the bath oscillators have been shifted by 



(19) 



still has a "quasiparticle" delta peak in the density of 
states 3m G^{uj) / tt, but of a reduced magnitude. It is su- 
perimposed on an "incoherent background" (see Fig. ||). 

For later use, we define the Fourier transform of the 
Green's function (without the constant energy shift): 



X exp 



iK{uj;K)= / die" 
'o 

dd, £M („ (1^1) + e{^)) (e-"^* - l) 



(20) 



For a bath with a spectrum that is stronger at low 
frequencies {a < 2), the Green's function may decay to 
zero for t —^ oo, which means that there is no delta peak 
any more in the density of states. 
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FIG. 2: The density of states for different bath spectra and 
temperatures. The quantity N{uj) = £smA"(a;)/7r, which cor- 
responds to the single-particle DOS obtained from the re- 
tarded Green's function, is plotted vs. frequency u/uic, for 
different exponents a of the bath spectrum at low frequencies 
{C{ui) ex. exp(— tj/ojc)) and for different temperatures. At 
T = 0, the DOS goes like uj°'~^ at low uj and vanishes for 
Lj < 0. In all cases displayed here, there is a delta peak at 
u — 0, except for a = 2 at T > 0. In the limit T — > oo the 
level shape becomes symmetric in all cases and the strength 
of the delta peak vanishes. 



The resulting single-particle density of states is pre- 
sented in Fig. H for different values of the exponent a 
characterizing the strength of the bath spectrum at low 
frequencies, both for zero and finite temperatures. 

Here, we will analyze in more detail the case of Nyquist 
noise, with a linear spectrum C(w) = Cqlo at small fre- 
quencies. 

If the temperature is finite, we may set n{uj) ~ T/lo for 
sufficiently low frequencies. The real part of the exponent 
(00) then becomes, in the long-time limit: 



27rCnK^T • t . 



(21) 



The prefactor in this expression defines the decay rate 
for the Green's function, which determines the width of 
the Lorentzian that arises in the density of states. 

At zero temperature the thermal excitation of the bath 
oscillators vanishes, so only the zero-point contribution 
remains. In that case, we have to specify the behaviour 
of the bath spectrum at high frequencies, because it be- 
comes important even at t — > oo: Choosing either a sharp 
cutoff C{uj) — Cqlu6{lUc — to) or an exponential decay 
C(w) = Cqw exp (— w/wc) results in the same long-time 
behaviour of the relevant integral: 
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FIG. 3: The density of states for Nyquist noise {a — 1) at 
T — 0, for different values of the power-law exponent 7. At 
finite temperature T, this goes over into a Lorentz peak of 
width (X T, see text. 



Here h{t) is a function that saturates to a constant 
value for i — > 00. The logarithm in the exponent leads to 
a power-law decay of the Green's function 



(23) 



with an exponent 7 proportional to the strength of the 
bath, the coupling and the momentum squared: 



(24) 



The detailed behaviour of the lineshape (not only the 
linewidth) then depends on 7, see Fig. Q. 

We will encounter such a power-law decay again in 
the discussion of dephasing for the Nyquist bath. It is 
similar to the power-law behay.iaur. found in the model 
of Quantum Brownian motiond'tI(E3. However, there are 
important differences between the Nyquist bath consid- 
ered here and the (stronger) Ohmic bath. These will be 
discussed in section IV D, 



IV. MANY-PARTICLE PROBLEM 

Now we consider the situation of many electrons on the 
ring, in the grand canonical ensemble at a given value of 
the chemical potential fi and at an arbitrary temperature 

T = f3-\ 



Grand canonical partition sum and persistent 
current 



doj 



C{lo) 



(cos [ujt) - 1) = -Co In [ujj^) + hit) . (22) 



As a preparation for the evaluation of various averages 
carried out below, we need the grand canonical partition 
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sum: 



/3H[fe}l 



(25) 



Here we have used the notation {pj} (which has been 
introduced above) for a given configuration of N parti- 
cles, trs in the second fine denotes the trace with respect 
to the bath oscillators. If we introduce an additional 
static magnetic flux <i>, this corresponds to a shift of the 
momenta pj ^ pj — e^/ Lc in all the formulas given be- 
low, such that the partition sum becomes dependent on 
$ as well. Note that we assume the effects of the charg- 
ing energy of the ring to be negligible here, unlike the 
treatment of cotunneling in subsection ^ 

The Hamiltonian can be split into a part con- 

taining only the momenta pj and another one which rep- 
resents bath oscillators which are shifted (depending on 
Pj) but still have the original frequencies (if we neglect 
the g'^cjP term like above): 



the kinetic energy T[{pj}] can become negative for large 
total momenta P if ^ > 1/iV. However, this only means 
that higher orders in ^ would have to be taken into ac- 
count (neglecting g^cf)"^ becomes invalid), as discussed al- 
ready in connection with the effective mass for a single 
particle (see Eq. [T^ ). 

Using this definition, we can rewrite Zgc, taking into 
account that the partition sum of the bath of shifted 
harmonic oscillators is equal to that of the unperturbed 
bath, Zho, a-nd therefore does not depend on {pj}'. 



Z,^ = ZhoY.^^^''T.' 



(29) 



N {p,} 



The partition sum without the contribution due to the 
harmonic oscillators will be denoted by Zgc from now on: 



^gc 



Zhc 



(30) 



We can simplify ( p9| ) further by rewriting Zgc in the fol- 
lowing way, using both the average (•)q and the partition 
sum with respect to the system of free electrons: 



N 

E 



p] - '^gp'f) 



H, 



bath 



2m 

T[{p,}]+H,ath[{Pj}], 



(26) 



where Hba 



th 



{pj}] is derived from the original bath 



Hamiltonian (G) by a shift in oscillator coordinates 



(27) 



and r[{pj}] is the residual "kinetic energy" term that 
only depends on the set of occupied momenta {pj}: 



(28) 



Here we have used the constant factor ^ defined in the 
discussion of the effective mass, see Eq. (p^). dif- 
fers from the single particle case by the appearance of the 
total momentum P = '^Pj, which makes it impossible to 
write this as a sum over the kinetic energies of individual 
electrons with renormalized masses. Physically, all the 
electrons are coupled to one and the same bath and this 
infiuences the center-of-mass motion. The dependence of 
T[{pj}] on the total momentum introduces a kind of sim- 
ple "effective interaction" between the electrons, which 
affects thermodynamic averages. For example, the aver- 
age occupation number for a given momentum does not 
follow the Fermi-Dirac distribution. Note that formally 



= (31) 

Here, we have kept only the lowest nonvanishing order 
in ^, since everything else would be inconsistent within 
the framework of our approximation. 

Now we can apply these results to the calculation of 
the persistent current of our system of electrons coupled 
to the fluctuating flux. 

The average current which flows for a given external 
magnetic flux $ is the derivative of the thermodynamic 
potential = —Tin Zgc with respect to the flux itself: 



= -{V[1-^N 



0^ 

(32) 



Here V — {P - 7Ve$/ {cL))/m is the total velocity 
operator for the electrons. We have replaced P in (|3ll ) by 
the expression valid in the presence of the external static 
flux. Obviously, the persistent current is reduced by a 
factor of about — The result does not depend 

on the details of the bath spectrum, only on f , which is 
well-behaved also for the case of the Nyquist bath. Recall 
that for values ^ > 1/iV the approximation used here is 
invalid, as explained above. 
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The magnitude of the reduction of the persistent cur- 
rent may be understood physically in the following way: 
If one imagines suddenly switching on the external mag- 
netic flux, an electric field pulse will be produced, which, 
at first, freely accelerates the electrons on the ring, lead- 
ing to a current which is proportional to the number N of 
electrons. This again produces a change in the magnetic 
flux which prompts a reaction of the bath (e.g. the exter- 
nal coil producing the Nyquist noise). The back-action 
onto the electrons deccelerates them, decreasing the ve- 
locity of each electron by an amount proportional to N 
and depending on the coupling strength between the ring 
and the external coil, which is contained in ^. This leads 
to the reduction factor 1 — obtained above. 

Thus, in our model, the reduction of the persistent cur- 
rent is similar in its origin to the appearance of an effec- 
tive mass. This can be seen most clearly by considering 
the special case of a "fast" bath, whose spectrum has a 
lower frequency cutoff or, at least, vanishes quickly with 
decreasing frequency. It follows the motion of the elec- 
trons adiabatically and is not able to lead to dephasing 
on long time scales (see the discussion of the two-particle 
Green's function below, as well as the cotunneling setup 
discussed in the last section). However, it still leads to a 
reduction of the persistent current, since the quantity ^ 
represents an integral over all frequencies. 



HO ■ e 



-PTl{p,}] 



(36) 



in order to define the average over the bath oscillators 
which are shifted depending on {pj}'- 



'l3H[{p,}] 



Zho^tb (e 



-/3Hbath[{p,}] 



(37) 



(See the discussion of Zgc in section [V A for a defini- 
tion of the quantities Zho and Hbath\{PiV\) 

With these definitions, expression (pq ) becomes: 



^Z{p^}{l-np[{pj}] 
{pj} 



Texp 



,^^t+t^ C dt'^{t') 

2m m Jo 



(38) 



In evaluating the average of the exponential, we need 
the expectation value of 4> which does not vanish in this 
case, since the oscillators are shifted by an amount pro- 
portional to the total momentum of the given configura- 
tion: 



B. Single-particle Green's function 

The single-particle Green's function is defined by 

iGp{t) = li'pt^o) m - (*po*Pt) ^("0 • (33) 

Evaluation of both expectation values proceeds in the 
same way, so we will only treat the first one here: 

^9c ^^""^^^ (e-'3^"e^^«**pe-^^"**t) . (34) 

N 

The trace is evaluated by summing over all configu- 
rations {pj} of N particles and using the fact that the 
bath does not change a given configuration. Therefore it 
is equal to: 

^ (1 - np[{pj}])trB (e-'3^[te}le*^[te}l*e-'^[{P^>'l*) . 
{pj} 

(35) 

Here {pj} is the configuration with one particle added 
in state p, and the prefactor is zero whenever that state 
is already occupied. Now we introduce the interaction 
picture with respect to Furthermore, we will 

use the partition sum for a given configuration. 



# = ..(^)^^^^=e(^E^^.j • (39) 

Apart from this, we can proceed exactly as before in 
the strictly single-particle case, see Eq. (^, in order to 
arrive at an exponent involving the thermal time-ordered 
correlator of for the unperturbed harmonic oscillators. 
Using this and (|3^) to evaluate (^), we have arrived at 
the desired result for one half of the Green's function, Eq. 
( |3^ ) . Proceeding analogously for the other half and using 
the definition ( po| ) for K{u;; k), the Fourier transform of 
the Green's function is given by: 

x{(i-n,[fe}])/^(.-|^4pP;f) 
-Mfe}]i^(-. + |l(i + O-V;f)}.(40) 

The sum runs over all configurations (of any particle 
number N) and P = ^Pj is the total momentum. 

From this result one can see that the lineshape of the 
DOS in the many-particle case is derived from the single- 
particle result K{uj; n). However, at finite temperatures. 
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it is the average over many such curves, each shifted by 
an amount cx pP that depends on the total momentum 
of the configuration. Still, the linewidth (or lineshape) 
does not depend in any essential way on the distance to 
the Fermi surface. 

We will now discuss the temperature-dependence of the 
linewidth that results from expression (pO|). The delta- 
peak which remains in K{u}) at w = for a "weak bath" 
(in this case: a > 2) is smeared over a certain range 
due to the average over configurations with different total 
momenta P. In lowest order with respect to ^, we may 
neglect the dependence of the probability distribution on 
the coupling to the bath. Then, the linewidth is obviously 
given by 



Slu — — p5P , 
m 



(41) 



where 5P — \J {P^)q is the spread in total momentum, 
calculated for the original free-electron system. We have 
{P'^}q = NmT and therefore a linewidth which increases 
with the square root of T: 



Suj cx ^Vt ■ 



(42) 



Note that the corresponding spread Sp — Slo/v in mo- 
mentum space is given by £,VNmT and can very well 
exceed the distance 2t:/L of the quantized momenta, in 
spite of the restriction <^ 1 (and also in spite of the re- 
striction \frrCr <^pf for the degenerate regime). There- 
fore, it is reasonable to speak of a linewidth, provided 
one does not resolve the quantized level-structure on the 
ring. 



time-arguments, but only for a perturbation acting at 
time and a density matrix evaluated at time t. Further- 
more, since the bath does not change the occupation of 
momentum states, the only nontrivial contribution arises 
from the following product of four electron operators, in 
which only two momenta p and p' appear, p refers to the 
hole that is created by the perturbation while p' belongs 
to the electron. 



-i-ffjvi,frt 



*T^r^,e-*^"'^';,,*p)(.43) 



N 



Inserting an appropriate basis of system states |{pj}) , 
using the interaction picture with respect to and 
carrying out the average of the exponential in the usual 
way, we arrive at the following result for the (half-sided) 
Fourier transform of Eq. (^) : 



^9^ X! "''P 



I ^ ( 2m ^ 2m ^ 



2m 



x"p[{pj}]{i - vKpj}]) 
— {p-p) + —{p -p)P\ 



) -(44) 



The notation is the same as for Eq. (^). The most 
important difference consists in replacing p hy p' ^ p in 
the factor k = gp/m which determines the strength of 
the decay. 

For the Nyquist case, we thus obtain a finite dephasing 
rate 



C. Dephasing: Two-particle Green's function 

While the decay of the single-particle Green's function 
in time is connected with every interaction process that 
changes the state of the electron or brings about random 
changes in its phase, it is not sufficient to know about 
this decay if one asks about dephasing. After all, there 
are situations where an electron interacts with a bath, 
such that its Green's function decays quickly but it is 
still able to show an interference pattern. This will hap- 
pen whenever the trace left by the particle in the bath 
is not enough to decide which path it Jaas gone, so that 
the different possibilities still interfereQ. Therefore, one 
must ask about the time-evolution of the density matrix 
(and, in particular, the decay of its off-diagonal elements) 
in order to study dephasing. Given a small initial per- 
turbation that creates a nonequilibrium situation, this 
time-evolution is determined by the two-particle Green's 
function in a linear-response calculation. 

The following calculation basically proceeds along the 
same lines as that given in the previous section, so we will 
keep it brief. For our purposes, we do not need the two- 
particle Green's function for arbitrary values of the four 



^1- = 2nTCo (^Yip'-pf , (45) 

which is proportional to the difference in momenta 
squared, the bath and coupling strengths and the tem- 
perature. The change in phase brought about by the fluc- 
tuating flux is proportional to p, so the phase-difference, 
whose variance appears in the exponent, goes like p' — p. 
To avoid confusion, we emphasize that there is no univer- 
sally applicable definition of a "dephasing rate" . In our 
case, we use this term to refer to the exponential decay 
of the two-particle Green's function as introduced above. 
Note that the dephasing rate can be small even if the 
decay rates associated with the (single-particle) Green's 
functions of the individual states are large. This is an 
example of the general behaviour mentioned above: It 
is reasonable that dephasing is strongest whenever the 
momenta of the two states, whose superposition is exam- 
ined, differ widely. Then, the bath, which couples to the 
momentum, can easily distinguish between these states 
even after a short time. 

The dephasing rate vanishes at T = 0. How- 
ever, even at T = the off-diagonal element of a density 
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matrix which initially describes a coherent superposition 
between momentum states p and p' decays to zero com- 
pletely in the limit t — s- oo. This decay proceeds with 
a power-law, as we have already observed for the single- 
particle Green's function, see Eq. (p3|). In the situation 
considered here, the exponent 7 is equal to 

l^(-)\p'-pfCo. (46) 
V m/ 

If the bath is sufRciently weak at low frequencies 
{C{uj) — Cqui" with a > 2), the decay of the off-diagonal 
elements in time saturates at a finite value, in contrast 
to the Nyquist case discussed above. Then the dephas- 
ing rate, defined as the prefactor of t in the exponential 
decay-law, is strictly zero, even at finite temperatures. 
This behaviour is related to the diagonal coupling be- 
tween system and bath. If a nondiagonal coupling were 
introduced, there would be transitions from excited elec- 
tronic states towards lower ones, accompanied by the 
spontaneous emission of a bath phonon. In that case, 
the decay rate in the single-particle Green's function and 
the dephasing rate would be nonzero also at T = but 
strongly dependent on the distance to the Fermi surface, 
due to the suppression of the density of final states for 
such transitions brought about by the Pauli-principle. 
Therefore, in a simple Golden rule calculation, these de- 
cay rates would vanish at T — when one approaches 
the Fermi surface. 

For a = 2, we find power-law dephasing only at finite 
temperatures, with an exponent proportional to T. 

D. Discussion of dephasing for the "Nyquist bath" 

The Nyquist bath is characterized by a fluctuation 
spectrum of flux and vector potential which is linear in 
uj (at zero temperature) , therefore leading to a spectrum 
for the electric field that behaves like uj^. This is ex- 
actly the spectrum of the zero-point fluctuations of the 
electric field in the vacuum. The main distinction be- 
tween those fluctuations and the Nyquist noise consid- 
ered here is that the latter leads to a force which is homo- 
geneous around the ring and therefore is compatible with 
the translational invariance of our one-dimensional sys- 
tem of electrons. Furthermore, its magnitude depends on 
the geometry and resistance of the external circuit pro- 
ducing the equilibrium current noise. Apart from these 
differences, we can use our understanding of the electro- 
magnetic vacuum fluctuations to discuss the effects of the 
Nyquist bath in a qualitative manner. 

In particular, free ballistic motion is not affected, since 
the radiation reaction force only acts on an accelerated 
charge. Therefore, the populations of the electronic mo- 
mentum eigenstates do not decay, as we have already ob- 
served. This is in contrast to the effect of a type of bath 
that leads to velocity-proportional friction, for example. 
The system is not ergodic, since the memory of the ini- 
tial conditions is not lost completely. In a basis other 



than the momentum basis, the off-diagonal elements of 
the density matrix show only partial decay. That there is 
some decay of the coherences is connected to the smear- 
ing of the position of the particle in the course of time. 
Usually, this effect is neglected in the discussion of dissi- 
pative quantum motion of a free particle under the influ- 
ence of a bath corresponding to electromagnetic vacuum 
fluctuations, since the ballistic expansion-pf an initially 
localized ensemble of particles dominateslla. In our case, 
it is important, since, for example, a superposition of two 
countcrpropagating plane waves on the ring will first form 
a standing wave pattern, whose visibility then gradually 
decreases. The fact that, at T = 0, the decay of the vis- 
ibility proceeds as a power-law can be understood most 
easily from the results of an oLipamiclassical analysis of 
the Lamb-shift due to WeltonEScJ: The vacuum fiuctu- 
ations of the electric field lead to a jitter of the electron 
position, such that the variance <^(5a;^) of its coordinate 
is given by the logarithm containing the ratio of an up- 
per cutoff-frequency (there taken to be the Compton fre- 
quency) and a lower cutoff (the characteristic frequency 
of electron motion around the nucleus). In our case, the 
lower cutoff frequency actually is given by the inverse of 
the observation time, such that (Jx^) oc In(wci)- For a 
superposition of plane waves of momenta ±k on the ring, 
the density matrix in position space contains cross-terms 
like exp(z2fca;), which, if averaged over Sx, give rise to a 
suppression factor exp (— 2fc^ {^^^)) ■ This leads to 

complete decay of the interference pattern even at T = 0. 
Note, however, that here we have been considering a su- 
perposition of excited states of the system and the decay 
of its coherences. In other problems of dephasing, such 
as those encountered in weak localization, one usually 
discusses the limit of zero-frequency response of the sys- 
tem to a small perturbation. A situation which comes 
closer to this kind of question will be discussed in the 
next section. 

The goal of the present work has been to analyze thor- 
oughly a model situation which shows some of the fea- 
tures important for dephasing at low temperatures, not 
to propose some experimental measurement setup. Still, 
we will now briefiy discuss the expected magnitude of the 
effect due to Nyquist noise in an external current coil. If 
the equilibrium current fluctuations are produced by an 
external coil whose circumference is similar to that of 
the Aharonov-Bohm ring (L) and which is placed about 
a distance L away, the dephasing rate is estimated to be 

\hcj \ c J R/Rk ^ ' 

Here Wi? is the Fermi velocity on the ring, R is the resis- 
tance of the external coil and Rk — h/e^ is the quantum 
of resistance. Both the square of the fine-structure con- 
stant in front of the expression and the ratio of the Fermi 
velocity to the speed of light render the effect very small 
under reasonable experimental conditions. 

Note that the fluctuations of the vacuum magnetic held 
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have a weaker power-spectrum, which leads to (AA)^ oc 
uj^ at T = 0, instead of {AA)^ oc w. The vacuum fluc- 
tuations of the electric field, however, do lead to a linear 
spectrum in the vector potential fluctuations (see discus- 
sion above). On the other hand, the electric field (at large 
wavelengths) is homogeneous only in free space, not with 
respect to its projection onto the ring, where it has a po- 
sition dependence cos {2t:x/L). Thus, the coupling is not 
diagonal in the momentum basis and is not included in 
our model. If one estimates the order of magnitude of the 
corresponding dcphasing rate (inelastic transition rate), 
one arrives at ~ {e^ /Tic) {vp / c)^ {vp / L). Simi- 

larly, one may estimate the strength of fluctuations due 
to shot-noise of the external current. In a situation where 
the external coil producing a static magnetic flux on the 
order of ^ $o is identical with that where the Nyquist 
noise originates, this leads to an effective dephasing rate 
°^ '^^{shot) ~ {e^ /Tic) {vp/c)vf/L, which may be much 
larger than that due to the Nyquist noise. Note that 
shot noise cannot be described by our model, since it 
is a nonequilibrium phenomenon and cannot be repre- 
sented by the usual bath of harmonic oscillators at low 
temperatures. However, the effects of shot noise would 
be reduced in a different geometry where a larger cur- 
rent (with correspondingly smaller relative magnitude of 
the shot noise) produces the same static magnetic flux 
through the Aharonov-Bohm ring. 



E. Relevance of the term quadratic in the flux 

In all of the preceding calculations, we have neglected 
the term g^cj)^ which appears in the kinetic energy of the 
particle but does not couple to the momentum. This ap- 
proximation has been necessary to use the well-known 
formula (exp {iX)) = exp (— {X'^) /2) for a Gaussian 
random variable X in the classical case or the analogous 
expression derived from Wick's theorem in the quantum 
case. (fP' is quadratic in the coordinates of the bath os- 
cillators, whereas </> is a linear bosonic variable (Gaus- 
sian random variable in the classical case). If the term 
g^cjp' /2m is kept in the Hamiltonian, the eigcnfrcquencies 
and normal coordinates of the bath of harmonic oscilla- 
tors are changed by an amount that depends on g but not 
on the particle momentum p. This becomes important at 
larger values of g, where the effective mass turns out to 
be m* — m{l + ^) instead of the value m* — m{l — ^)^^ 
derived without the 0^-term, see Eq. (17). Note that 
this difference persists also in the thermodynamic limit 



oo. However, the qualitative behaviour of the 



bath spectrum at low frequencies is not changed for the 
spectra C{u)) oc uj°' with a > 1 that have been consid- 
ered here. Therefore, this term is unimportant for the 
qualitative conclusions about dephasing, although it can 
change quantitative results for larger g and does change 
the single-particle Green's function. The latter involves 
a change in particle number, so that an additional 4>^ is 




FIG. 4: The tunneling setup discussed in the text. 



introduced into the Hamiltonian. 



V. AHARONOV-BOHM INTERFERENCE IN 
COTUNNELING THROUGH THE RING 

In the following, we will discuss the influence of the 
bath on the Aharonov-Bohm (AB) effect, i.e. on the flux- 
dependence of the transport current through the ring. 
We consider tunneling into and out of the ring, taking 
place at two electrodes to the left and right of the ring 
(see Fig. ^). A tunnehng situation is the appropriate 
one for our model, since attaching current leads would 
severely alter the system. Note that independent (se- 
quential) tunneling only probes the density of states at 
the two contacts. Therefore, in order to observe the AB 
interference effect, we have to consider a Coulomb block- 
ade situation, in which any electron tunneling into (or out 
of) the ring will enhance the total energy by the charg- 
ing energy of the ring which is much larger than the bias 
voltage V and the temperature T. In such a case, traus.- 
port through the ring is possible only via cotunnelin^il, 
i.e. a two-step process involving a virtual intermediate 
state belonging to a different number of electrons on the 
ring. A strong dependence of the tunneling current on 
the external magnetic flux, with a complete suppression 
at $o/2 due to destructive interference, is visible only 
in the "elastic cotunneling" contribution, where the elec- 
tronic state of the ring is left unaltered in the process. It 
is linear in the bias voltage and will dominate the inelas- 
tic contribution at low temperatures and for small bias 
voltages (see the discussion at the end of this section and 
Ref. |l|). 

Now consider cotunneling at T = under the influ- 
ence of the bath. The semiclassical analysis of Aharonov- 
Bohm interference given in a preceding section is not ap- 
plicable for V — > 0, since it assumes the electron can 
emit or absorb an arbitrary amount of energy. In the 
quantum-mechanical calculation, suppression of interfer- 
ence is due to the electron leaving a trace in the bath that 
permits, at least in principle, to decide which of the two 
arms of the ring the electron has traveled. This involves 
a transfer of energy between electron and bath. The bath 
spectrum determines the amount of bath oscillators able 
to absorb the small energy < eV which can be emitted 
by the electron. Therefore, it is to be expected that de- 
phasing at zero temperature is suppressed for — > due 
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to the energy conservation constraint. This will be con- 
firmed by the calculation described in the following, al- 
though there are renormalization effects that change the 
strength of the tunneling current away from the point of 
perfect destructive interference, $ — ^o/2. 

The tunneling process starts from a situation in which 
the ring is occupied by the equilibrium number of elec- 
trons (depending on the value of a gate voltage) and the 
Fermi seas in the left and right electrode are filled up 
to Fermi energies that differ by the bias voltage, eV, see 
Fig. ^. Throughout the following discussion, we will as- 
sume the Fermi energy in the left electrode to be the 
larger one of the two (and eV > 0). In the final state, an 
electron has appeared above the right Fermi sea, leaving 
behind a hole in the left electrode. Although we want to 
consider the situation where the electronic state of the 
ring has not changed in the end, the final state of the 
bath may be different. The intermediate state is char- 
acterized by an extra electron (or extra hole) present on 
the ring and some arbitrary state of the bath. Using 
standard second-order Fermi's Golden rule, the tunnel- 
ing rate is obtained by summing over all intermediate 
states (dividing by the proper energy denominator) and 
all final states whose energy equals the initial energy: 



2-E 

/ 



E,, — Ei 



6{Ef-E,). (48) 



Here = T'"-\-T^ is the sum of the tunneling Hamil- 
tonians belonging to the left and the right junction, while 
the energies and eigenstates refer to the unperturbed 
Hamiltonian that includes everything besides tunneling. 
In particular, it includes the coupling between electrons 
and the bath, as well as the kinetic energies of electrons 
in the electrodes. 

At this point we would like to emphasize that using 
Fermi's Golden rule for the calculation of the cotunncl- 
ing current does not in itself mean taking into account 
the interaction between bath and system only in a per- 
turbative way. The intermediate states being summed 
over in Eq. (Q) are the exact eigenstates of the full sys- 
tem of electrons on the ring coupled to the fluctuating 
flux. In this sense, the coherence properties of the ring 
as a whole (including the bath) are tested by the cotun- 
neling process. Applying the Golden rule in this context 
is roughly comparable to using the Kubo formula in a 
linear-response calculation of, e. g., the weak- localization 
magnetoconductance, which does not automatically im- 
ply a perturbative description of the dephasing processes 
either. Still, there is an important difference: In our con- 
text, we essentially deal with a scattering situation, such 
that energy conservation holds at least for the initial and 
final states of the complete process. This will be seen to 
be important for our conclusions about the strength (or 
absence of) dephasing. What is neglected in our calcu- 
lation are any equilibrium correlations between the state 
of the electrons in the leads and the state of the ring as 




FIG. 5: Energy diagram for cotunneling through the AB-ring, 
at $ = $o/2. The initial, final and two possible intermediate 
states are indicated (see main text). The charging energy has 
to be added to the single-particle energy of the intermediate 
state shown here. 



a whole (including the fiuctuating fiux). 

Before performing the calculation in the presence of the 
bath, we will briefiy describe how the destructive inter- 
ference at $o/2 appears in this formula, in the situation 
without fiuctuating fiux. In such a case, the intermedi- 
ate state V refers solely to the electronic state k on the 
ring, which is occupied by the additional electron in the 
course of tunneling. The final state / is determined both 
by the state A, which is unoccupied in the left electrode 
after the tunneling process, and the state A, where the 
electron ends up in the right electrode (see Fig. |5|). 

For simplicity, we will assume tunneling to take place 
only between two points, for example from a point at 
the tip of the left electrode to an adjacent point xl on 
the ring. 



(49) 



and likewise for the right electrode, is a complex- 
valued tunneling amplitude. Such a description will be 
appropriate as long as the extent of the relevant region in 
which tunneling can take place is less than a wavelength. 

The sum over intermediate electronic states k on the 
ring then contains the following contribution which de- 
scribes an electron going onto the ring from the left elec- 
trode and leaving through the right electrode: 



ek+ Ec - e\ 



(50) 



Here the sum over k is to be taken only over unoccu- 
pied single-electron states on the ring. '5 refers to single- 
electron wavefunctions on the ring and on the electrodes. 
Apart from the contribution listed here, there is another, 
completely analogous, contribution which belongs to the 
situation with an extra hole on the ring in the interme- 
diate stateEJ. Note that for the purposes of our discus- 
sion we will not distinguish between the charging energies 
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belonging to the electron- and hole-processes (assuming 
them to be of about the same magnitude). 

Perfect destructive AB-interference at an external 
static magnetic flux of $ = $o/2 arises only for an even 
number of electrons on the ring. In this case, the energies 



C(W) 



_ 1 / 27r $ 



(51) 



of the unoccupied states are pairwise degenerate, for 
k+ EE n27r/L and fc_ = (1 - n) 2n/L, see Fig. |. There- 
fore, the energy denominators for k+ and k- are the 
same, while the wavefunctions in the numerators produce 
a phase shift of exp {i {k+ — k-) L/2) = —1 between the 
two possibilities, leading to complete cancellation of all 
terms in the sum. The same applies to the sum over oc- 
cupied states (for the situation with an extra hole in the 
intermediate state). 

Dephasing will, in general, "wash out" this perfect de- 
structive interference. After taking the modulus squared 
of the sum of amplitudes given above in Eq. (|50|), which 
we briefly denote by A{k), one obtains "classical" prob- 
abilities like |A(A:_|_)| but also cross-terms of the form 
A* {k+)A(k-). A bath coupling to the electronic motion 
will affect these terms differently, if it is able to "distin- 
guish" between the momenta and k_ . Usually, the 
cross-terms are suppressed. Then, the different contri- 
butions cannot cancel any more. Note that, away from 
perfect destructive interference, we have to expect an in- 
fluence of the bath on the magnitude of the tunneling cur- 
rent under any circumstances, since mere renormalization 
effects like a change in the effective mass of the electrons 
will be important. This is why we concentrate on the 
special case of $o/2- Even in that situation, the interfer- 
ence minimum could vanish in a rather trivial way due 
to renormalization effects, if one chose a bath coupling 
asymmetrically to the two arms of the ring (or rather to 
left- and right-going momenta), thereby leading to dif- 
ferent transmission probabilities. This is not the case in 
our model. 

At this point, we can give a simple counting argument 
to arrive quickly at the voltage-dependence of the "inco- 
herent" contribution to the cotunneling rate F, which is 
produced by the fluctuating flux and leads to a nonvan- 
ishing current at $o/2. The sum over initial electronic 
states on the left electrode is carried out over a region of 
extent eV. The probability of emission of a bath phonon 
is proportional to the bath spectrum oc w", and we have 
to integrate this from to the maximum energy of the 
electron, which is again of order eV. This yields a volt- 
age dependence F cx 1/"+^ of the incoherent contribution 
to the tunneling current (see Figs. ^ and H). However, 
we want to analyze this in the following using a different, 
nonperturbative scheme, thereby making contact to the 
Feynman- Vernon influence functional formalisnJj, which 
is the "workhorse" of many dephasing calculations. 

To evaluate Eq. ( |4^ ) in the presence of the bath, 
we rewrite the sum over intermediate states as a time- 
integral: 




FIG. 6: Cotunneling in the presence of the fluctuating flux: 
Emission of a bath "phonon" of frequency uj leads to an in- 
coherent contribution to the tunneling current, where the de- 
structive interference between the two paths shown here is 
lost. The probability of emission depends on the bath spec- 
trum C{ij). 



ttT ttT „oo 

T.ir^-' dt{f\H^e-^^"'^^^'H^\^) (52) 
- hi Jq 

k 

In the second line, which replaces Eq. (|5^) in the pres- 
ence of the bath, we have split off the contribution due 
to the electronic states. This has been possible because 
only the tunneling operators can change the electronic 
state, while the bath couples diagonally to the electronic 
momentum eigenstates. Furthermore, we have confined 
ourselves to the process with an extra electron in the in- 
termediate state, a completely analogous contribution for 
an extra hole has to be added. H[k] is the Hamiltonian 
for a given configuration consisting of an extra occupied 
state k over the original Fermi sea on the ring. It only 
acts on the bath Hilbert space, where 0^ refers to the 
ground state of the bath prior to the tunneling event and 
is an arbitrary final state which the bath goes into 
after the cotunneling process is finished. In our nota- 
tion, the sum of electronic kinetic energies is included in 
H[k] as well, whereas the charging energy Ec has been 
taken into account separately. The matrix elements of 
the tunneling Hamiltonians T^^^ are taken between the 
electronic states A, k and A (compare Eqs. (^ ) and (po|)). 
Eq is the ground-state energy of the ring, including the 
bath. 

After taking the modulus squared of the sum given 
above, we arrive at the following contribution to the tun- 
neling rate F at zero temperature: 



X 



2^ E iT-^k>Tt>x)iT-^u<Tt<xr (53) 

A,A,fe>,fc< 











X J2(^<ir<)\nSiEfn - i?o - (ex - 6,))(/^|x>(r>)) , 
where fc^^*"^ denote unoccupied states on the ring. 
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There are three analogous contributions besides the one 
shown here, in which the tunneling takes place in a dif- 
ferent order (e.g. the process may start by an electron 
tunneling out of the ring, leaving a hole behind, etc.). 
Note that a similar expression arises in the derivation 
of the "P(i?)" -thftory of a tunnel junction coupled to a 
dissipative bathE3E3. 

In the preceding formula, the last line can be 
viewed as a kind of "generalized influence functionaF 
FIt-' ,t'^,uj = e\ — ex]. It is equal to the overlap between 
bath states x^''^' which have been time-evolved out of 
0^ under the action of for some time t-'^-^K 

In contrast to the usual influence functional, the time 
of evolution may be different for the two states and the 
overlap is taken only with respect to bath states at an 
excitation energy uj (which must equal the energy emit- 
ted by the electron). This difference is due to the fact, 
that in our problem the energy conservation constraint 
must be taken care of, since the electron cannot transfer 
an arbitrary amount of energy to the bath. This clearly 
shows why a single-particle calculation using the usual 
influence functional must fail when the amount of energy 
available is limited due to low temperatures or low bias 
voltages. This problem has also been discussed in Ref. ||, 
where the authors have used physical arguments to drop 
certain "zero-point" contributions to the dephasing rate 
obtained in a single-particle calculation. The normal in- 
fluence functional is recovered by integrating over all pos- 
sible energy transfers and setting r> — t"^: 



F[r]^{x^ (r) \x> (r)) 



+ 00 



dujF [r, r, u] . (54) 



If no bath is present or its spectrum has a lower cutoff 
which is larger than the energy available to the electron, 
energy conservation leads to a; = 0, such that the final 
and initial bath states coincide: = 0^. Then, F is 
a product of a factor depending only on fc> and another 
one, depending only on k^. In this case, the sums over 
/j>(<) niay be carried out separately, like before, and the 
terms will cancel again (for $o/2), provided the bath 
couples equally to k+ and fc_ (see discussion above) . Al- 
though there is definitely no dephasing in this case, the 
magnitude of the tunneling current may be changed for 

7^ $0/2, due to the afore- mentioned renormalization 
effects. 

The Fourier transform (in lo) of the generalized influ- 
ence functional may be written as follows: 



F[r>,r<,r] 



27r 



(x<(r<)|e-^«^«lnx>(T>)) . (55) 



It can be repr esented as a Keldysh time-ordered expec- 
tation valueESt3, apart from a prefactor exp(— ii?o(T> — 
r<)): 




L AB ring R 



FIG. 7: Schematic "space-time" diagram showing the Keldysh 
contour which runs from to r and back again. The interac- 
tion operator which couples to the bath is nonvanishing only 
when there is an extra electron on the ring (thick lines), either 
in state or fc^ (see main text). Tunneling processes are 
indicated. 



Here, Vi{s) = gk(j){s)/m + k^/{2m), with k = k — 
27r$/(<I>oi), where <I> is the additional static flux. Vj cou- 
ples the additional electron in state k to the bath (and 
also incorporates the kinetic energy). We have k = k^ 
if s is on the forward time-branch and < s < r^, 
while A; = fc^, if s is on the backward time-branch and 
T -I- — < s < T + . For all other times, V/ van- 
ishes. This is represented graphically in Fig. 0. Note that 
Vi is taken in the interaction picture with respect to the 
bath coupled to the original Fermi sea on the ring. Once 
more, we have neglected the term which turns out 
to be unimportant for the bath spectra considered below. 

Using Keldysh time-ordering and a linked cluster ex- 
pansion (Wick's theorem), we can represent Eq. ( ^6| ) as 
an exponential containing double time-integrals involv- 
ing the Kcldysh-time ordered correlation function of the 
bath operator (j): 



exp 



i j^ds (Vi(s)) 



(57) 



X exp 



K 



K 



dsi <b ds2(TKSViisi)SVi{s2) 



(56) 



We have set SVj = Vi — {^ij = gkcjy/m. The princi- 
pal steps involved in the evaluation of this expression are 
demonstrated in Appendix where we show how the 
usual Caldeira-Leggett influence functional for a bath of 
harmonic oscillators can be derived very efficiently using 
this method. The resulting exponential couples the mo- 
menta fc^ and fc^ and may therefore lead to dephasing. 

From now on, we again consider baths which are char- 
acterized by a power-law spectrum at low frequencies, 
C{uj) cx uj" with an exponent a > 1. Remember that the 
case a = 1 represents fluctuations of the magnetic flux 
produced by Nyquist noise of an external current loop. 
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FIG. 8: Schematic behavior of the cotunnehng rate as a func- 
tion of static magnetic flux and bias voltage. In the ideal 
situation, the rate vanishes at $0/2 (dotted line), while it 
rises as a power of the bias voltage due to the incoherent 
contribution resulting from the fluctuations of the flux (thick 
line). 



For these bath spectra, it is sufficient to carry out an 
expansion of the cotunneling rate to leading order in the 
coupling strength g. The part of the resulting expression 
which couples fc^ and A:^ is seen to lead, after summation 
over all electronic states A:^, fc^. A, A, to an "incoherent" 
contribution that washes out destructive interference but 
which is suppressed for low bias voltages, as expected. 

At zero temperature, the ratio of this incoherent cur- 
rent at $ = $0/2 to the normal elastic cotunneling cur- 
rent which flows at $ = is found to be given by the 
following approximate expression (up to a constant of 
order 1): 



2 2 

9 



eV 



(1 



eV 



)C{uj)duj 



1 

6? 



(58) 



The expression inside the brackets can be interpreted 
as the variance of the fluctuating energy of a single- 
particle level on the ring. However, it is to be evalu- 
ated taking into account only the fluctuations up to the 
frequency corresponding to the bias voltage and using a 
weight factor I—lu/ eV which favors low-energy transfers 
uj. We have already pointed out that the cutoff at eV 
is a simple consequence of energy conservation. For a 
power-law bath spectrum C(w) cx w", the integral yields 
a voltage dependence cx V"'^^, so the incoherent tunnel- 
ing current goes like Note that Se = hvp/L refers 
to the single-particle level-spacing on the ring. The qual- 
itative behaviour of the cotunneling rate as a function 
of both external static flux and bias voltage is shown 
in Fig. (^). (We remark again that the incoherent cur- 
rent due to external Nyquist noise would be too small to 
be actually measurable under reasonable experimental 
conditions, compare the discussion at the end of section 
iVD| ). 

Thus, we see that suppression of destructive interfer- 



ence does not show up in the linear conductance. In 
this sense, the fluctuations do not lead to dephasing in 
the linear transport regime at zero temperature. For the 
Nyquist case a = 1, the exponent of V is the same as 
that for inelastic electronic cotunneling processes in a sks.- 
tem with a continuum of intermediate electronic statesciJ. 
Bath spectra with a > 1 obviously lead to an even weaker 
decrease in the visibility of the interference minimum at 
low bias voltages. Note that formally inserting a = — 1 
in r cx I/"+2 would lead to an incoherent contribution 
to the linear conductance even at T = 0. However, this 
case is not of interest here, since it cannot be produced 
by a fluctuating magnetic flux and it is not covered by 
the approximations made in our calculation (in particu- 
lar dropping the 0^-tcrm). It would correspond to the 
strong force fluctuations of an Ohmic Caldeira-Leggett 
bath used in the description of quantum Brownian mo- 
tion. 

For bias voltages eV smaller than the single-particle 
energy spacing Se on the ring, dephasing is merely due 
to the coupling to the fluctuating flux. At higher volt- 
ages, the inelastic cotunneling processes become impor- 
tant. In these, one electron tunnels into the ring, while 
another electron goes out at the opposite electrode, thus 
leaving behind a particle-hole excitation on the rin^il. 
Since all the corresponding final states are different, their 
contributions to the cotunneling current sum up incoher- 
ently. Therefore, like dephasing produced by the bath, 
they also lead to a nonvanishing contribution to the tun- 
neling current at $o/2, where, ideally, one should have 
perfect destructive interference. The number of possibil- 
ities to create a particle-hole excitation with an energy 
of at most eV is cx (eV/Se)"^, if we assume eV ^ Se. In 
that regime, the ratio of the incoherent current contribu- 
tion due to the external bath to the electronic inelastic 
contribution is given by the bracket in (^), multiplied 
by {Se/ (EceV))'^ . The electronic inelastic contribution 
will be the dominant one. 

Finally, let us discuss finite temperatures. Without the 
bath and as long as T <C Se, only the Fermi distributions 
in the electrodes get smeared, which does not affect the 
tunneling current, if one takes into account that now the 
tunneling processes do not only lead to an electron trans- 
port from left to right but in the other direction as well. 
The presence of the bath will introduce some tempera- 
ture dependence for the incoherent current contribution 
in this regime, since at finite temperatures the tunneling 
electron can not only emit an energy quantum into the 
bath but also absorb a thermal bath excitation. There- 
fore, the energy to transferred to the bath now can be 
negative as well. There is no restriction on the amount 
of energy an electron can absorb, so there is no cutoff eV 
for negative to. At positive energy transfers, the proba- 
bility of spontaneous emission into the bath (cx C(w) in 
Eq. (p8|)) now has to be multiplied by n{uj) + 1, where 
n{Lo) is the Bose distribution function (induced emission). 
At negative w, this is replaced by n(|a;|), since only ab- 
sorption of thermal excitations (not of vacuum fluctua- 
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tions) is possible. 

Taking into account the thermal smearing of the 
electrode Fermi distributions, the balance of left and 
right-going tunneling currents and the induced emis- 
sion/absorption mentioned just before, we have to per- 
form the following replacements in (p8|): 



r du:{l-^)C{u:) 



+ OC 



duj W{p, w, eV) {n{\uj\) + e{uj)) C{\lu\) . (59) 



Note that the factor multiplying W in the integral cor- 
responds directly to the function P{E) which occurs, in 
the theory of tunneling in a dissipative environmentE3c3. 
The function W itself represents an integral over the av- 
erage of energies in the left and right electrodes, e = 
(ea + Ca) /2; at a fixed energy transfer uj = tx — e^- 



Wi(3,iu,eV) 



ieVy' J defL{l-fR)-fR{l-fL)e- 



(60) 

Here h = fp (ea - eV/2) and fn = fp (e^ + eV/2) are 
the Fermi distributions in the two electrodes. The factor 
of exp is due to the fact that for an electron going 

from right to left the energy transferred to the bath is 
—ui, so the ratio between the values of 't.(I'jjI) + [uj) 
at positive and negative frequencies appears, which is 
just this factor. Note that, in the low-temperature limit 
(/? oo), the function W becomes 1 — uj/eV for uj < eV 
and vanishes for lo > eV. This reproduces the left-hand 



side of Eq. ( |59| ) as a special case. 

Using this formula, the incoherent tunneling current is 
found to be enhanced by a temperature-dependent con- 
tribution (X T"+iy. 

At T > 5e, one would have to take into account the 
thermal averaging over different electronic configurations 
on the ring (still at a fixed particle number determined by 
charging energy and gate voltage). The perfect destruc- 
tive interference at $o/2 depends on the presence of an 
electronic configuration which is symmetric in the occu- 
pancy of equal-energy states having fc > and k < (see 
discussion above). The thermal average includes other 
configurations as well and therefore leads to a suppres- 
sion of the destructive interference in the elastic tun- 
neling current, even without the bath. Furthermore, the 
electronic inelastic contribution is also enhancecLat finite 
temperatures and becomes linear in the voltagfd. 



between findings in different physical situations is not 
straightforward, since there is no general definition of a 
"dephasing rate". E.g., in recent work on an Aharonov- 
Bohm ring containing a quantum dot capacitively cou- 
pled to the Ohmic fluctuations in a gatet2l, the authors 
found that the coupling suppresses the magnitude of the 
persistent current flowing in the ring and induces fluctu- 
ations of this current. The effect persists even at T = 
and was interpreted as dephasing at zero temperature. 
Formally, the coupling assumed in the setup of Ref. |l^ 
is nondiagonal in the system's eigenbasis, whereas it is 
diagonal in our case. This may indeed lead to a weaker 
tendency towards dephasingin our model (see the discus- 
sion at the end of Section y) . Still, we also find both a 



reduction of the persistent current (subsection |IV A| ) and 
fiuctuations (caused by the vector potential fluctuations, 
since the momentum is conserved), regardless of the de- 
tails of the bath spectrum. Nevertheless, this apparently 
does not imply dephasing in our interference setup (co- 
tunneling transport situation), which has been our main 
concern here and which has no analogy in the work of 
Ref. ^ The transport effect depends in an essential 
way on the long-time dynamics of the system and, there- 
fore, on the low-frequency behaviour of the bath spec- 
trum (similarly to the Caldeira-Leggett model or some 
aspects of the spin-boson model), while the reduction 
of the persistent current does not. Therefore, we have 
not interpreted the reduction of the persistent current 
in terms of zero-temperature dephasing, although it cer- 
tainly constitutes a suppression of an interference-related 
phenomenon. We believe that the suppression of interfer- 
ence in a cotunneling setup considered here may be more 
directly related to mesoscopic transport interference ex- 
periments. 

Our model is certainly quite removed from the discus- 
sions about low-temperature dephasing in weak localiza- 
tion. There, dephasing by electron-electron interaction 
in an extended disordered system weakens the coher- 
ent backscattering effect, which is an interference phe- 
nomenon that is robust against both thermal and impu- 
rity averaging, in contrast to the destructive interference 
considered in our model. Since the electrons inside the 
metal interact all the time, it is not obvious whether that 
situation is in any way analogous to the kind of "scatter- 
ing" situation considered here for the cotunneling trans- 
port. 



VII. CONCLUSIONS 



VI. DISCUSSION 

Since we do not find dephasing in our model at zero 
temperature (under the circumstances specified above), 
the reader may wonder how our result is related to some 
apparently contradicting conclusions in the literature. 
First of all, we would like to point out that a comparison 



We have analyzed a simple model of a fluctuating mag- 
netic flux threading an Aharonov-Bohm ring and dis- 
cussed its effects on equilibrium properties, such as per- 
sistent current and tunneling density of states, as well 
as transport properties such as the two-particle Green's 
function and the cotunneling current through the ring. 
Particular emphasis has been put on the question of 
dephasing and low-temperature behaviour. There are 
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important qualitative differences depending on the low- 
frequency behaviour of the bath power spectrum. For 
the case of Nyquist noise and an arbitrary initial su- 
perposition of momentum states, an exponential decay 
of off-diagonal elements of the density-matrix in time at 
T > goes over into "power-law dephasing" at zero tem- 
perature. However, if one probes the coherence proper- 
ties of the electronic motion on the ring by checking for 
the possible suppression of destructive AB-interference 
in a cotunneling setup, no such suppression is found in 
the linear transport regime at T = 0. This is because 
the possibility for the electron to leave a trace in the 
bath is diminished due to the energy conservation con- 
straint. Our calculation serves as an illustrative example 
of the difference between the "optics" type of interference 
experiments (employing single particles), in which the 
semiclassical approximation and/or the usual Feynman- 
Vernon influence functional may be applied to discuss 
dephasing, and the linear-transport interference experi- 
ments encountered in mesoscopic physics, in which spe- 
cial care has to be taken in the analysis of dephasing at 
low temperatures. 
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average over xo has to be performed at finite tempera- 
tures: 



J{x<,x>) = {U{t,0\x<)xo \U{t,0\x>)xo) 



Xo 



(Al) 



Here the time-evolution operators depend on the sys- 
tem trajectory x{-) via the interaction term V{t) = x{t)4> 
in the Hamiltonian. Using the interaction picture with 
respect to the bath Hamiltonian Hb, we can explicitly 
write down J in the following form, with (anti-)time- 

ordering symbols T (T): 



J = Zg^tr 



(A2) 

Here Zb = tr exp{—(3HB)- Now we can use Keldysh 
time-ordering to abbreviate this formally: 



J = ZgHr 



(A3) 



Here xk{s) is equal to x^{s) (or a;<(s)) if s lies on 
the forward (or backward) part of the Keldysh contour 
that runs from to i and back again. The advantage 
of this formal rearrangement is that the application of 
Wick's theorem (leading to a linked-cluster expansion) 
now becomes very simple. Wc immediately obtain for J: 



APPENDIX A: DERIVATION OF 
CALDEIRA-LEGGETT INFLUENCE 
FUNCTIONAL USING KELDYSH 
TIME-ORDERING AND WICK'S THEOREM 

The influence functional for a system (variable x) cou- 
pled linearly to a linear bath (0) is usually derived in 
the pathpigtegral picture, by "integrating out" the bath 
variablestHd'cl. This can be done because the bath consists 
of a set of uncoupled harmonic oscillators. However, the 
calculation is usually quite cumbersome, although the re- 
sult is simple enough and involves only the real and imag- 
inary parts of the bath correlator (4>{ti)<f>(t2)^ ■ Here 

we present a derivation based on Keldysh time-ordering 
and Wick's theorem. The oscillators of the bath and 
their action for an external driving force never have to 
be considered and therefore this is probably the shortest 
route to the well-known Caldeira-Leggett influence func- 
tional. Modifications of this approach are applicable in 
more complicated situations as well (compare the main 
text). 

The Feynman- Vernon influence functional is the over- 
lap between the two bath states which result from the 
action of two different (fixed) system trajectories x''{-) 
and x^{-) onto the same initial bath state xo- A thermal 



exp 



K 



K 



dti <p dt2 (TK(j){tl)(bit2)) XKitl)xKit2) 



, (A4) 

The brackets (•) denote the thermal average with re- 
spect to Hb ■ Now we can translate back the exponent by 
taking into account all four possible combinations of the 
two times on the forward or backward time paths: (To 
keep the notation short, we use subscripts for the time 
arguments) 



dU 



2 } a;.^ -f 



H(p2 ) xfx^ 



5201 ; X 



.>x<}. 



(A5) 



This can be simplified further by noting that the real 
part of all the four different correlators appearing here is 
the same, since it is symmetric in the time arguments: 

23?e^(^i(^2^ = (^^4>i,4>2^^- Therefore, the real part of 
the exponent is given by: 



dti / dt2^e 



(A6) 
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It defines the imaginary part of the influence action 
Sinfi[x'^ jX"^] and describes dephasing and heating due 
to the fluctuations of the bath variable 4>. 

Treating the imaginary part is only slightly more dif- 
ficult. We have 



This defines the real part of the infiuence action and 
describes friction and renormalization effects (e.g. ef- 
fective mass). Not e that we can bring the real part of 
the exponent, (A6), to a similar form by cutting off the 
t2-integral at ti and dropping the factor 1/2 in front of 
the expression, since there the integrand is symmetric in 

^1,^2- 



(A7) 



3to(T0i02) = sgn{ti - t2)'^m (4'i(j)2 



3m ( T(j}i(j)-2 



In order to get rid of the sgn(ti ~ 12), we split the 
double time-integral into one part where t2 < ti and 
one with t2 > ti. In the latter part, we interchange 
integration variables ti and t2- Then we obtain for the 
imaginary part of the exponent: 



dti 



dt2 5m (0102 ) {xi - xf) {x2 + X2) 



(A8) 



In this way we have arrived at the well-known influence 
functional for a system coupled linearly to a linear bath. 
As usual, a couphng of the form f{x)(j) just leads to a 
replacement x 1— > f{x) in the final expression, and some- 
thing like fj{x)4>j leads to a sum of the corresponding 
influence actions, if the (h^ are uncorrelated. 
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